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Infrared behavior of the Landau gauge gluon propagators is studied numerically in the 3d SU(2) 
gauge theory on the lattice. A special accent is made on the study of Gribov copy effect. For this 
study we employ an efficient gauge-fixing algorithm and a large number of gauge copies (up to 280 
copies per configuration). It is shown that, in the deep infrared region, the Gribov copy effects are 
very significant. Also we show that, in the infinite-volume limit, the zero-momentum value of the 
propagator does not vanish. 
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I. INTRODUCTION 

The nonperturbative (first principle) numerical 
computation of the field propagators is important for 
various reasons. There are scenarios of confinement 
based on infrared behavior of the gauge dependent 
propagators. In particular, in the Gribov-Zwanziger 
(GZ) confinement scenario [l|, the Landau-gauge 
gluon propagator D(p) at infinite volume is expected 
to vanish in the infrared (IR) limit p 0. On the 
other hand, a refined Gribov-Zwanziger (RGZ) sce- 
nario proposed recently allows a finite nonzero 
value of -D(O) . Another reason is that the nonpertuba- 
tive lattice calculations are necessary to check the re- 
sults obtained by analytical methods, e.g., the Dyson- 
Schwinger equations (DSE) approach which uses trun- 
cations of the infinite set of equations. The DSE scal- 
ing solution predicts that the propagator tends to zero 
in the zero-momentum limit [6j, [7[ in accordance with 
the RGZ-scenario. At the same time, the DSE decou- 
pling solution allows a finite nonzero value of 
D(0) in conformity with RGZ-scenario. 

The 3d SU(2) theory can serve as a useful test- 
ground to verify these predictions in a simplified, in 
comparison with the Ad case, setting. Furthermore, 
3d theory is of interest for the studies of the high- 
temperature limit of the 4c? theory. 

The most theoretically attractive definition of the 
Landau gauge is to choose for every gauge orbit a 
representative from the fundamental modular region 
[12j . i.e. the absolute maximum of the gauge fixing 
functional F{U) (see the definition in Section ILT1L 

The arguments in favor of this choice are the follow- 
ing : a) a consistent non-perturbative gauge fixing 
procedure proposed by Parrinello-Jona-Lasinio and 



Zwanziger (PJLZ-approach) [TH, [l4| presumes that 
the choice of a unique representative of the gauge orbit 
should be through the global extremum of the chosen 
gauge fixing functional; b) in the case of pure gauge 
U(l) theory in the weak coupling (Coulomb) phase 
some of the gauge copies produce a photon propagator 
with a decay behavior inconsistent with the expected 
zero mass behavior [l5|-[l7||. The choice of the global 
extremum permits to obtain the physical - massless - 
photon propagator. It should be noted that, for prac- 
tical purposes, it is sufficient to approach the global 
maximum close enough so that the systematic errors 
due to nonideal gauge fixing (because of, e.g., Gribov 
copy effects) are of the same magnitude as statistical 
errors. We follow here this strategy which has been 
checked already in many papers on 4d theory studies 
for both SU{2) [H-[22[ and SU(3) [H, gauge 
groups. 

The three-dimensional SU{2) theory has been re- 
cently studied in [26 31]. The evidence has been pre- 
sented that the propagator has a maximum at mo- 
menta about 350 MeV and that -D(O) does not con- 
verge to zero in the infinite- volume limit. The prob- 
lem of Gribov-copy effects was addressed in [3(| • The 
Gribov-copy effects for the gluon propagator were 
found for small momenta. We will show in this pa- 
per that these effects were underestimated. 

As compared to the previous version of this article, 
we take into account the lattices 80 3 and 96 3 and ob- 
tain better statistics on lattices of smaller size. For 
this reason, we make the respective improvements in 
our plots, and draw more definite conclusions. 

In Section |TT] we introduce the quantities to be com- 
puted. In Section [HI] some details of our simulations 
are given. In Section IIVI we discuss the effect of im- 
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proved gauge fixing and present our numerical results. 
Conclusions are drawn in Section fVl 



and i 
one obtains 



—gaA/2. For zero- momentum propagator 



II. THE GLUON PROPAGATOR: 
DEFINITIONS 

We consider cubic Lx Lx L lattice A with spacing 
a. To generate Monte Carlo ensembles of thermalizcd 
configurations we use the standard Wilson action 
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(1) 

where (3 = A/g 2 a, fx is a vector of length a along the 
fj,th coordinate axis and g denotes dimcnsionful bare 
coupling. U xfl € SU(2) are the link variables which 
transform under local gauge transformations g x as fol- 
lows: 



Ux» 4 U°„ = gtU x tf x+lk , g x e SU{2) . (2) 
We study the gluon propagator 

D%{d) = ^ E exp(iqx){Al(x + y)Al(y)}, (3) 
where the vector potentials are defined as follows [321 ]: 

^ = E<T = i(^-^)' ( 4 ) 

6=1 Z 9a 

and the momenta take the values = 2irn fl /aL, 
where runs over integers in the range — L/2 < < 
L/2. The gluon propagator can be represented in the 
form 



9L 3 



(7) 



H=l 6=1 



In the weak-coupling infinite-volume limit, D(p) = 

—z. In this work, we deal with the quantity (referred 
pi 

to as the gluon propagator) 1 



D{p) = D(p)/p; 



(8) 



In lattice gauge theory the usual choice of the Lan- 
dau gauge condition is [32[ 

(dA)(x) = 1 E (A^x) - A„(x - A)) = , (9) 

which is equivalent to finding a local extremum of the 
gauge functional 



F u(9)= ^E^ Tr "it 



3L 3 ^ 2 



(10) 



with respect to gauge transformations g x . The 
manifold consisting of Gribov copies providing local 
maxima of the functional (|10[) and a semi-positive 
Faddeev-Popov operator is termed the Gribov region 
f2, while that of the global maxima is termed the fun- 
damental modular region (FMR) rcfl. Our gauge- 
fixing procedure is aimed to approach F. 



III. DETAILS OF THE SIMULATION 
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where p^ = — sin and p 2 — 53^=1 v\- When 
p ^ 0, we arrive at 
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where 



(5) 



° b M) = E cos ( 9a; ) u U x ^ ( 6 ) 

xEA 



x£A 



In this work, we are going to demonstrate that the 
gluon propagator in the deep infrared region can be re- 
liably evaluated only when the effects of Gribov copies 
are properly taken into account. We make simula- 
tions at a fixed gauge coupling on various lattices us- 
ing gauge-fixing algorithm which was already succes- 
fully employed in the 4<i theory at both zero [2(| Hl[ 
and nonzero [22l . |25| temperature. There are three 
main ingredients of this algorithm: powerful simulated 
annealing algorithm, which proved to be efficient in 
solving various optimization problems; the flip trans- 
formation of gauge fields, which was used to decrease 
both the Gribov-copy and finite-volume effects 20 
|22| ; simulation of a large number of gauge copies for 



1 Note that the function D(p) (not D(p)\) is normalized similar 
to the analogous function used in |2S| and [3CH . 
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each flip sector in order to further decrease the effects 
of Gribov copies. 

We perform Monte Carlo (MC) simulations at (3 = 
4.24 for various lattice sizes L. Consecutive configura- 
tions (considered to be statistically independent) were 
separated by 200 sweeps, each sweep consisting of one 
local heatbath update followed by 15 microcanonical 
updates. In Table U we provide the full information 
about the field ensembles used in this investigation. 



L 


fl>meas 


f^copy 


ai[fm] 


Prmn [GeV] 


32 


800 


96 


5.38 


0.230 


36 


1037 


160 


6.05 


0.204 


40 


1032 


160 


6.73 


0.184 


44 


717 


160 


7.39 


0.167 


48 


1425 


160 


8.08 


0.153 


52 


1085 


160 


8.74 


0.141 


56 


796 


160 


9.43 


0.131 


64 


709 


160 


10.8 


0.115 


72 


910 


280 


12.1 


0.102 


80 


557 


160 


13.5 


0.092 


96 


438 


280 


16.1 


0.077 



TABLE I: Values of lattice size, L, number of measure- 
ments n me as and number of gauge copies n copy used 
throughout this paper. Spacing is a — 0.168 fm. 

The features of the gauge-fixing methods [lj| used 
in our study are as follows. Firstly, we extend the 
gauge group by the transformations (also referred to 
as Zi flips) defined as follows: 



u x . 



if 



and 



a, 



Ux.p, otherwise 



which are the generators of the Z\ group leaving the 
action ([TJ invariant. 

Such flips are equivalent to nonperiodic gauge trans- 
formations. A Polyakov loop directed along the trans- 
formed links and averaged over the 2-dimensional 
plane changes its sign. Therefore, the flip operations 
combine the 2 3 distinct gauge orbits (or Polyakov loop 
sectors) of strictly periodic gauge transformations into 
one larger gauge orbit. 

The second feature is making use of the simu- 
lated annealing (SA), which has been found compu- 
tationally more efficient than the use of the standard 
overrelaxation (OR) only [H |H, Hi]. The SA al- 
gorithm generates gauge transformations g(x) by 
MC iterations with a statistical weight proportional 
to exp(3V Fu\g]/T) . The "temperature" T is 
an auxiliary parameter which is gradually decreased 
in order to maximize the gauge functional Fjj[g] . 
In the beginning, T has to be chosen sufficiently 
large in order to allow traversing the configuration 



space of g{x) fields in large steps. As in Ref. [Jfj . 
we choose Ti n i t = 1.5. After each quasi-equilibrium 
sweep, including both heatbath and microcanonical 
updates, T is decreased with equal step size. The 
final SA temperature is fixed such that during the con- 
secutively applied OR algorithm the violation of the 
transversality condition 



3 

max \ 

x,a — ^ 



(u a Jx)-u a Jx~f,)) 



< e 



(11) 



decreases in a more or less monotonous manner for 
the majority of gauge fixing trials until the condition 
(fTTj) becomes satisfied with e = 10~ 7 . A monotonous 
OR behavior is reasonably satisfied for a final lower 
SA temperature value Tfi na i = 0.01 [33j. The number 
of temperature steps is set equal to 3000. The final- 
izing OR algorithm using the standard Los Alamos 
type overrelaxation with the parameter value ui = 1.7 
requires typically a number of iterations of the order 
O(10 3 ). 

We then take the best copy out of many gauge fixed 
copies obtained for the given gauge field configuration, 
i.e., a copy with the maximal value of the lattice gauge 
fixing functional F as a best estimator of the global 
extremum of this functional. 

In what follows, we consider three gauge-fixing 
methods. 

The first employs both the SA-OR algorithm and 
the Zi flips. Using the SA-OR algorithm, we generate 
8n c Gribov copies (n c copies in each Z\ sector) and 
find the copy giving the maximum of the functional 
(|TU|) . This copy is referred to as the best copy ("6c ") 
and it should be mentioned that we also find "the best 
sector" for each starting configuration. The version of 
the Landau gauge obtained by this method is labelled 
FSA ("Flipped Simulated Annealing"). 

The second method employs the SA-OR algorithm, 
whereas the Zi flips are not taken into consideration. 
In this case, we choose "the best copy" from n c config- 
urations, that corresponds to a random choice of the 
Z\ sector. We name it "SA gauge-fixing" (it is anal- 
ogous in some sense to the "absolute gauge fixing" in 

In order to demonstrate the effect of Gribov copies, 
we also consider the gauge obtained by a random 
choice of a copy within the first Gribov horizon (la- 
belled as "/c " — first copy). That is, in the third case, 
we take the first copy obtained by the SA algorithm 
and do not take care of fundamental modular domain 
at all. This version of the Landau gauge is analogous 
in some sense to the "minimal Landau gauge" in [30|). 

Information on our simulation procedure is also 
shown in Table Q] The scale is set by the relation 
a^/d « 0.376 for /3 = 4.24 (see formula (65) in (H), 
where a w (440 MeV) 2 is the string tension (which is 
assumed the same as in the four dimensions). 
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FIG. 1: D(0) as a function of n copy . The meaning of 
'copy-sector' and 'sector-copy' is explained in the text. 



IV. GRIBOV COPY EFFECTS AND LARGE 
L BEHAVIOR OF D(0) 

To demonstrate the Gribov copy effects, we show 
in Fig. Q] the dependence of D(0) on the number of 
gauge copies n copy for L — 36 lattice. As one can 
see, the Gribov copy influence is very strong (at least, 
for p = 0), and by no means can be considered as a 
" Gribov noise" . 

For the upper curve in this Figure ('copy-sector') 
first 20 gauge copies belong to the first flip-sector, 
next 20 gauge copies belong to the second flip-sector, 
etc.. Evidently, the dependence on the number of 
copies belonging to the same flip-sector is rather weak 
(apart from the first one), and the main Gribov copy 
effect comes from different flip-sectors. This is demon- 
strated also by the lower curve ('sector-copy') where 
we have used another enumeration of gauge copies, 



''copy 



1 corresponds to the 1st copy of the 1st 
flip-sector, n copy — 2 means that first copies of the 
first and second flip-sectors are considered, etc.. With 
increasing volume n copv — dependence is changing as 
we will explain below in discussion of Fig. [3l 

To compare Gribov copy effects for different values 
of L and p ^ 0, we define the Gribov copy sensitivity 
parameter A(p) = A(p; L) as a normalized difference 
of the fc and be gluon propagators 



A( P ) 



D /C (p) - D bc (p) 
D bc (p) 



(12) 



where the numerator is the average of the differences 
between fc and be propagators calculated for every 
configuration and normalized with the be (averaged) 
propagator. 



In Table [TT| we show the values of the parameter 
A(p) for different values of L and four momenta : p = 
0, Pmin, V^Prnin and V^Pmin- It is interesting to note 
that values of this parameter do not demonstrate the 
tendency to decreasing with increasing lattice size L 
for all momenta under consideration. In particular, at 
zero momentum the Gribov copy effect is estimated to 
be between appr. 25% and 35% (taking into account 
the error bars) for all values of L. 



L 


p = 


Pmin 


pmin 


pmin 


32 


0.318(14) 


0.107(7) 


0.043(5) 


0.014(6) 


36 


0.329(13) 


0.107(6) 


0.039(4) 


0.013(6) 


40 


0.289(13) 


0.104(7) 


0.039(5) 


0.031(6) 


44 


0.308(16) 


0.100(8) 


0.042(5) 


0.026(7) 


48 


0.251(11) 


0.100(5) 


0.046(4) 


0.029(5) 


52 


0.285(13) 


0.091(6) 


0.043(4) 


0.031(6) 


56 


0.273(15) 


0.116(8) 


0.042(5) 


0.023(7) 


64 


0.294(17) 


0.100(8) 


0.049(6) 


0.027(8) 


72 


0.261(16) 


0.121(7) 


0.052(5) 


0.044(6) 


80 


0.232(20) 


0.102(10) 


0.062(7) 


0.031(8) 


96 


0.218(21) 


0.115(12) 


0.074(8) 


0.055(9) 



TABLE II: Values of A(p) for the FSA Landau gauge prop- 
agators at p — 0, p = Pmin, P = V2p m in and p = V3p m in- 

Therewith, for a given value of L, the parameter 
A(p) decreases quickly with an increase of the momen- 
tum. This observation is in agreement with the obser- 
vations made earlier for the four-dimensional SU(2) 
theory dJ. 

We have attempted to estimate the infinite- volume 
limit of the zero-momentum propagator D(0; L), i.e., 
the limit L — > oo. 

One can apply various fit-formulas for this purpose, 
e.g., C1+C2/L; C2/L Ca ; ci + C2/i Ca , etc., and many 
of them fit nicely if the values of L are not very large 
(as in our case). However, calculations [28| on the 
lattice with L — 320 and j3 = 3.0 have shown that the 
first fit-formula is supposed to be the preferable one 
(at least, in the minimal Landau gauge). Therefore, 
following |2S| we apply the fit-formula 



D(0;L) =ci +c 2 /L 



(13) 



to determine the L — > 00 limit of D(0;L). 

In Fig. [2] we show our values of D(0; L) calculated 
for a) fc (which is the same for SA and FSA meth- 
ods); b) be SA method (i.e., without flips) and c) 
be FSA method. Broken lines represent fits accord- 
ing to Eq. (|T3l) . 

We confirm that in the L — > 00 limit the value of 
D(0) differs from zero. This is in agreement with the 
statement made in [28| and is not in conformity with 
M and EL 
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Fig. [2] shows another interesting phenomenon : the 
Gribov copy influence survives even in the thermo- 
dynamic limit L — > oo. Indeed, the infinite volume- 
extrapolation of D^ c (0) differs from infinite- volume 
extrapolation of D bc (0) (both for SA and FSA gauge 
fixing algorithms). 

To illustrate this phenomenon in a more explicit 
way, we calculated also the averaged (over all configu- 
rations) difference between fc and be propagators nor- 
malized to D bc (p; L — oo) (in what follows it will be 
refered to as W(p)). In Fig.[3]we show the dependence 
of W(Q) on the inverse lattice size both for SA and 
FSA algorithms. One can see that both algorithms 
predict nonzero difference between fc and be values 
of the propagators, this difference being rather big 
(~ 15%) even in the thermodynamic limit. Note 
that both SA and FSA algorithms give the coincident 
(within errorbars) results in the L — > oo limit. 

For small volumes W(0) for SA algorithm, Wsa(0), 
is close to zero while Wfsa{0) is at its maximum. 
This corresponds to strong effects of flip-sectors seen 
in Fig. [T] With increasing volume Wsa(0) is increas- 
ing indicating the increasing role of the copies within 
given flip-sector. In opposite, decreasing of Wfsa(0) 
with increasing volume implies that the role of the 
flip-sectors reduces. In the infinite volume limit the 
two curves in Fig. [3] converge. This means that in 
this limit one randomly chosen flip-sector is sufficient 
or, in other words, all flip sectors are equivalent. Note 
that on our largest volume the effect of the flip-sectors 
is still dominating over the effect of the copies within 
one sector. 

In Fig. fj] we compare the L-dependence of W(0) 
shown in the previous Figure with that for two 
nonzero values of momenta: \p\ — 150 MeV and 
\p\ = 250 MeV (all calculated for FSA). The respective 
values of the propagators for fixed physical momenta 
were obtained by interpolation of the original data to 
necessary value of momentum. As it was expected, 
Gribov copy effects decrease quickly with increasing 
\p\. However, they are expected to be not small in the 
deep infrared region, i.e., in the \p\ — > limit. This 
can be essential for the calculation of, say, screening 
masses in 4d theory at nonzero temperature where the 
IR-behavior of the gluon propagator is important. 

In Fig. [5] we compare the momentum dependence 
of the propagators D(p) calculated on 36 3 , 72 3 and 
96 3 lattices. Qualitatively its momentum dependence 
agrees with that obtained earlier [3(3, HH Ht| . But in 
contrast to these papers the finite-volume effects are 
very small and can safely be neglected at momenta 
p > 0.4 GeV. Note that the propagator has a maxi- 
mum at nonzero value of momentum \p\ ~ 400 MeV. 
Therefore, the behavior of D(p) in the deep infrared 
region is inconsistent with a simple pole-type depen- 
dence. 
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FIG. 2: D(0) as a function of L for fc ( the same for SA 
and FSA), be SA (without flips) and be FSA. 
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FIG. 3: W(0) as a function of L for FSA and SA. 



V. CONCLUSIONS 

In this work we investigated numerically the Landau 
gauge gluon propagator D (p) in the three-dimensional 
pure gauge SU(2) lattice theory. The main goal was to 
study the approach of the zero- momentum propagator 
D(0) to the thermodynamic limit L — > oo. We have 
employed eleven lattice volumes from L = 32 to L = 
96. All calculations have been made at /3 = 4.24 (a = 
0.168 fm). 

Special attention in this study has been paid to the 
dependence on the choice of Gribov copies. To this 
purpose we have generated up to 280 gauge copies 
for every configuration. Our be FSA method provides 
systematically higher values of the gauge fixing func- 
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FIG. 4: W(p) as a function of L for FSA for three different 
values of \p\ . 
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FIG. 5: The momentum dependence of D(p) for three 
lattices. 



tional as compared to the fc FSA and be SA methods. 
The main results are the following. 



1. In the limit L — > oo the value of D(0; L) differs 
from zero (in agreement with [28j). 



2. The Gribov copy influence is very strong in the 
deep infrared region (especially for D(0)), and 
by no means can be considered as a "Gribov 
noise" . 

Moreover, our analysis shows that the Gribov 
copy effects remain substantial (up to ~ 15% for 
zero-momentum propagator) even in the ther- 
modynamic limit L — ¥ oo. The choice of the 
efficient gauge fixing procedure is of crucial im- 
portance in this study. 



3. The finite- volume effects for the gluon propaga- 
tor at momenta greater than 400 MeV are neg- 
ligibly small. 
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